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We show that quintessence, when it is described by a tachyonic field, can amplify a tiny primordial 
gradient generating a preferred direction in the sky. In its simplest realization, this mechanism only 
affects the Cosmic Microwave Background fluctuations at the quadrupole level. We briefly discuss 
how higher multipoles can also be affected, once the full structure of the quintessence potential is 
taken into account. 
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I. INTRODUCTION 

The formulation of the cosmological principle pQ co- 
incides with the birth of modern scientific cosmology. 
Over the last decades, homogeneity and isotropy of the 
Universe have been tested to increasing degrees of ac- 
curacy. In particular, the radiation of the Cosmic Mi- 
crowave Background (CMB) is one of the most sensitive 
probes to test the isotropy of the Universe. The observed 
temperature of the CMB is observed to be uniform at ze- 
roth order (once the dipole is subtracted): this result 
led to the formulation of the theory of primordial infla- 
tion. Moreover, approximate statistical isotropy appears 
to hold even for the small fluctuations of the CMB tem- 
perature. 

In spite of these striking results, several analyses of the 
CMB fluctuation maps, starting from [2113] (see, e.g., [3] 
for a more recent analysis), have shown the existence of 
anomalies associated to some degree of breaking of sta- 
tistical isotropy. Even though a clear consensus on the 
subject is still absent {e.g., see jS] for an up to date dis- 
cussion) intensive work has been done in the past also 
as a response to the theoretical challenge of naturally 
generating a preferred directions in the sky. The mod- 
els proposed usually rely either on early-Universe or on 
late-Universe mechanisms. While the former lead to an 
intrinsic anisotropy in the primordial spectrum of met- 
ric perturbations, the latter produce anisotropy in the 
observed CMB via anisotropic Sachs- Wolfe effect. Early- 
Universe mechanisms include Hj (see however [7]) and [5], 
that rely on vector fields or [9], that rely on spinor fields, 
as well as [TQl [TT], that rely on primordial gradients. 
Late-Universe mechanisms can invoke magnetic fields [H] 
or anisotropics in the dark energy equation of state [T51 - 
[15] . Despite those many attempts, it is fair to say that 
breaking statistical isotropy usually requires strong and 
aesthetically unappealing assumptions. 

In the present work we show a (relatively) natural 
mechanism that can give rise to a preferred direction in 
the CMB sky. This mechanism is a hybrid of early and a 
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late Universe ones and does not rely on vector fields or on 
unusual properties of the dark energy sector, but rather 
on a simple model of scalar dark energy, characterized by 
a non-trivial, tachyonic (V"((f) < 0) potential. 

Models of tachyonic quintessence frequently appear 
in the literature. For instance, if the quintessential 
scalar is described by a pseudo-Nambu-Goldstone Boson 
(pNGB) [TS], then half of the extrema of the potential 
V(<j>) oc 1 + cos(0//) are tachyonic. Another scenario of 
tachyonic quintessence was proposed in [T7], where the 
magnitude of the tachyonic mass is related to the height 
of the potential at its maximum. 

Our main assumption is that a tachyonic field <j> has 
a small primordial gradient. As the Hubble parameter 
drops below y/\V"(jtj>)\, the initial gradient is amplified 
by the "pull" of its tachyonic potential, effectively con- 
verting a small primordial isocurvature perturbation into 
a larger curvature mode. In first approximation, this 
leads to an "ellipsoidal" Universe [T2] which could ex- 
plain the low quadrupole CMB amplitude observed in 
both COBE and WMAP data [H], if the preferred di- 
rection in the quintessence correlated with that of the 
primordial quadrupole. 

As we will see, our mechanism will be at work for a 
tachyonic mass a few times larger than A/Mpj. This is, 
in particular, the case for a pNGB with value of the axion 
constant / slightly sub-Planckian. 

The primordial gradient in <fr can be a remnant of the 
epoch that preceded the last bout of inflation, provided 
that such last bout were sufficiently short (this argument 
will be discussed in greater detail in section [v]). In 
this respect, our scenario could shed some light on 
pre-inflationary dynamics, similarly to the situations 
studied, e.g., in [TH] and especially [TO]. In our case, a 
constant gradient can be the dominant remnant of a more 
general inhomogeneous initial value of the tachyonic 
field, i.e., we can expand <fi w k + k% z p + n 2 z p + . . . - 
for sake of simplicity, we will only consider powers of the 
z coordinate, though the argument for the dominance 
of the linear term would still hold, were more generic 
terms considered - at Ni = N ohs + SN e-foldings before 
the end of inflation (N i, s corresponds to the time when 
the current cosmological scales left the horizon and z p to 
the physical distance along the z direction). We assume 
that each term Kk z p had given an equal contribution 
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to the energy in (f> at that moment. Today, the scales, 
that left the horizon N ^ s + 5N e-foldings before the 
end of inflation, are still outside the horizon by a factor 
of e SN , i.e., they correspond to physical distances of 
the order of e SN Hq 1 . Assuming no (other) significant 
evolution in tj> [20] , the assumption that each term Kk z p 
be of the same order implies that Kk ~ e~ h5N Hq . As 
a consequence, inside the horizon z p ss Hq 1 , each term 
Kk Zp contributes like e" kSN to the energy in 4>. This 
implies that the terms with the lowest powers in z p will 
give the largest contributions to the metric perturbations 
at sub-horizon scales. For this reason, our analysis will 
be focused on the approximation where, in the early 
Universe, <f> = K\z with (f>(z = 0) being set equal to zero 
by an appropriate choice of the origin of the z coordinate. 

The paper is organized as follows. In section [TT] we 
set up our system and solve the corresponding Einstein 
equations. Sections [TTT] and |TV| contain the effect of our 
anisotropic metric on the redshift and in particular on 
the CMB quadrupole anisotropy. Section [V] contains a 
discussion of the magnitude of the primordial gradient 
responsible for the anisotropy. Finally, in the concluding 
section we will discuss how our mechanism could also 
lead to alignments of the higher multipoles. 



II. THE METRIC 

We consider the cylindrically symmetric cosmological 
metric 

ds 2 = -dt 2 + a(t, zf (da; 2 + dy 2 ) + b{t, zf dz 2 . (1) 
The non- vanishing components of the Einstein tensor are 



G f± = 



2 a' 


b' 


1 


b 2 ~a~ 


b 




Gyy = 


-a 2 


(- 


a' 6 




of 




- 2 


•t 


a b 




a 



a 2 a b 
a z a b 

1 a' b' la" ab d V 



b 2 a b b 2 a ab a b 



G Z z — b ( „ „ „ 2 . 

b A Gr a z a I 



(2) 



where a dot denotes a derivative with respect to t and a 
prime a derivative with respect to z. 

We assume the presence of three contributions to the 
stress-energy tensor: pressure-less dust with energy den- 
sity p, a cosmological constant A and a scalar field <f> with 
canonically normalized kinetic term and potential V{4>) 
- in principle one can always consider A to be a part of 
V(4>) however we prefer to keep these two components 
separate as we will keep separate the cosmological con- 
stant problem and the breaking of statistical isotropy. 



The components of the stress energy tensor of <j> read 



T X x — ^yy ® 



2 V 



1 

262 



^ - Vtf) - A 



T, 



1 

262 



1 " ^' 2 - V(cp) - A 



(3) 



and throughout the paper we will consider (j) to be a 
tachyon with potential V(cf>) = —m 2 4> 2 /2. The stress en- 
ergy tensor for dust has the general form T^ v — pu^Un, 
with u 1 " V^ L u v = 0. 

Let us now solve the above equations by assuming 
that the non-isotropic, non-homogeneous scalar field is 
the source of a perturbation around the Fricdmann- 
Robertson- Walker metric of a ACDM Universe. As a 
consequence, we write a(t, z) = ao(t) [1 + Sa(t, z)] and 
b(t, z) = a Q (t) [1 + 5b(t, z)], where 5a < 1 and <56 < 1; 
the Friedmann equation for the background = has 
exact solution 
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where Hq is the current value of the Hubble parameter 
and Q\ ~ 0.3 from observation. 

We can assume that the dust velocity field is equal 
to (1 + 5u°, 0, 0, 8u z ) with Su° and 8u z being first order 
quantities. Thus the geodesic equation for reduces 
to SvP = and [a^5u z y — 0. Choosing Su z = as 
t —> 0, we obtain, up to second order in the perturbations, 
?i M = (1, 0, 0, 0). Then conservation of the stress-energy 
tensor for dust implies that p = po/(o 2 6), where po is 
the dust density at a fixed time. 

Turning on the scalar field <fi, we assume that, as de- 
scribed in the introduction (see also [HJ HU the 
lowest terms in the expansion in powers of the coordi- 
nate z give the largest contributions to the anisotropy. 
Then, without loss of generality, we can set 4>(z = Q) = 
(remember that <f> = corresponds to the maximum of 
the tachyonic potential), so that 



<P(t, z) = K X {t)z + 0{z 2 ) . 

At leading order, the Klein-Gordon equation 
(j)" ' /oq — m 2 = decomposes as 
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We then expand 5a(t, z) — 6an(t) + 8a,2(t) z 2 + 0(z 3 ) 
and Sb(t, z) = Sb (t) +<56 2 (£) z 2 + 0(z 3 ) (it is easy to see 
that, since at this level of approximation <j) is odd in z, 
da and 5b must be even in z). The (tz) Einstein equation 
- at first order —2 6a' = <p' <p/Mp^ - gives rise to 
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The (zz) Einstein equation reads 
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The only non-trivial part of this equation is the one in- 
dependent on z, as the rest is redundant with respect to 
the two previous equations: 



Lastly, we compute the correction to the scale factor 
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along the z direction, 8bo, from (11) 
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The (it) Einstein equation reads, at first order in 8a 
and (56 and using the background equations, 
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where the only term under control in the z n expansion is 
the one in z° (since 8a is determined with 0(z 3 ) accuracy, 
the presence of the term 8a" constrains the accuracy of 
this equation to be 0(z) hence, at this level of approxi- 
mation, 8b can be determined only in its z° component). 
Therefore the only equation we obtain from the (tt) com- 
ponent is 
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where Sip® = 2 Sao + £60 ■ 



The explicit expressions of Ki(t), 8ao(t) and 8bo(t) are 
the following. The solution of eq. using the expres- 
sion Q for ao(t), reads 
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where we defined the dimensionless quantities r = 
^y/TTfiHo t and /1 2 = ^m 2 /{^n^H^), and where the inte- 
gration constant R\ is defined so that Ki(0) = R\ Hq Mp\. 

It is straightforward to solve for 8a 2 from eq. ^ using 
the found solution for cf>(t, z) at the desired order: 
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where we have imposed 8a 2 {t — > 0) = 0. 

Sa Q can be computed by integrating eq. d9j): 
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III. THE PERTURBED REDSHIFT 

An observable quantity related to the metric found in 
the previous section is the angular dependence of the 
CMB radiation induced by the integrated Sachs- Wolfe 
effect on our system. We compute it in this and in the 
following section. Photons coming from the last scatter- 
ing surface move along null geodesies: 



d£+ 
dA~ 







(16) 



where A is the affine parameter along the worldline 2+ (A) 
of the photon, k^ — da+/dA denotes the tangent (null) 
vector to the photon worldline. 

It is convenient to introduce the conformal time 77, so 
that the metric reads 

ds 2 =al{ri) [-dr, 2 + (1 + 2 8a(r], z)) (dx 2 + dy 2 ) + 

+(1 + 2 56(77, ^))dz 2 ] . (17) 

This metric will allow us to factor out the dependence 
of the redshift £ [53] on the scale factor so to effectively 
simplify the calculation. Let us perform the conformal 
transformation 

ds 2 = a 2 (?7) ds 2 , dA = a 2 (r]) dA , k" = a^y 2 ^ . 

(18) 

It is easy to show by taking advantage of fe M being a null 
vector, that k^ satisfies 
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where T^ p are the Christoffels symbols of the metric ds 2 , 
which is a first order perturbation around the Minkowski 
metric. 

The photon worldline, 2+ (A), can be expanded in a 
perturbation series 

a+(A) = x^CXj+Sx^iX) , fc"(A) = fc£(A)+<5fc"(A) , (20) 

where the subscript B denotes a background quantity 
and 8 denotes the first-order perturbation on it. 

Since the background component of ds 2 describes a 
Minkowski spacetime, then the background null vector 
fcn is constant. By appropriately scaling A we impose 



k% = -1. We have k% 



x£'(A) = (- 



1, n), where n is 
a unit 3-vector, as we will eventually choose A = for 
the observer of the CMB photon, with A increasing as we 
go back in time. Given the symmetries of our spacetime, 
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we can locate the observer at (xb)o — {vb)o — 0, while 
(zb)o = zq remains arbitrary, so we have 

Vb(X) = t] -X , (21) 
x B (X) =n x X, y B (X) =n v X, z B (X) = z + n z X . 

The redshift of a source with four- velocity Ug, as seen 
by an observer whose four- velocity is u , is defined as 
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It is convenient to introduce a "conformal" redshift (, 
defined by adding a bar to all the quantities in eq. ( 22 ) 
and defining = 1 (jj)u^: 
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where we have used the fact that ao(?7o) = 1- We assume 
that the spatial components of Ug and Uq can be treated 
as first order quantities, in which case C is also a first 
order quantity. By using Uq s — (1, vo.s), where vo is 
the peculiar velocity of the observer and vs of the source, 
we compute C to be 



C ~ (vs - Vo) 



5k (X s ) + Sk°(X = 0) 



(24) 



where we have used the fact that the photon was emit- 
ted at A = As and is observed at A = 0. The quantity 
5k°(Xs) — Sk°(0) is obtained by integrating the geodesic 
equation for k° 
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where the first term corresponds to the intrinsic motions 
of the source and of the observer. The second term is 
associated to the integrated Sachs- Wolfe effect on our 
spacetime; the next section will be dedicated to a thor- 
ough discussion of this term. 



IV. CMB ANISOTROPIES 

The effect of our anisotropic space on the CMB spec- 
trum is given by the integrated Sachs- Wolfe (iSW) effect: 
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Using (21 1, we trade A for rj, moreover, since Sa and 5b 
rapidly vanish as rj — > 0, we can extend the integration 
to rj = (i.e., As = rjo). We remind that Sa ~ Sao(rj) + 



60,2(7)) z 2 and 8b m Sbo(rj). Since 60,2 and <5&2 appear on 
the same footing in the equations above, and since <5&2 is 
undetermined at this level of approximation, we neglect 
also S a 2 in what follows. Thus we have 
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where n z = cosi9. 

The above expression can be integrated explicitly 

CiSW « <5a (?7o) + (56 (t7o) - SaoiVo)) cos 2 , (29) 

in which we have used 8ao(rj = 0) = 6bo(i] = 0) = 0. It is 
assumed that the direction of alignment of the gradient 
coincides precisely with the z axis. In order to find the 
equivalent expression when the gradient of <p is directed 
along a direction (8, $), we rotate the coordinate system 
(x, y, z) by an angle $ around the z axis and by an angle 
around the x axis. We thus obtain 

Cisw(0, <?)» 5a (rj a ) + (8b (r) ) - Sa (rio)) x (30) 
x [sin sin cos(<I> + </?) + cos $ cos 0] 

We follow, e.g., [H] to find the effect on CMB. 
If T* is the CMB temperature at decoupling, then 
the observed temperature will be T($, p) = T*/(l + 
C(0, ip)) ~ T, acMB (1 - C(^) <P))> where a C MB is the 
scale factor at decoupling. Then the average observed 
temperature is (T) = J dip d cos $ T '(1? , ip)/ATr and the 
anisotropy 5T{&, ip)/T = 1 - T(0, <p)/f = C iSW (0, <p) - 
J dip d cos CiSW (0, <p)/4nc. Its decomposition in spheri- 
cal harmonics is given by 
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An explicit calculation allows to show that the only 
non-vanishing contributions from the anisotropic inte- 
grated Sachs- Wolfe effect are 



4 S 2 W = \/^ (Sb (rj ) - Sa (m)) sin 2 0e 2j * , (32) 



iSW 
l 21 



ASW 
l 20 



(Sb Q (Tjo) - 5a (r, )) sin 20e* $ , 



3 J 5 (SboM - Sa (r] )) (1 + 3 cos 20) 



15 

1 Ptt 



with 02,-1 = — a 21 an d ^2,-2 = ^22' 

These contributions have to be summed to the intrin- 
sic quadrupole momentum of the CMB. It is possible to 
see that our contribution can help explaining the low 
quadrupole observed in CMB data, since it effectively 
gives rise to an "ellipsoidal" Universe |12) . 

Let us quickly review the argument of [T2]: we de- 
note by a\ m (with a\ _ m = (— l) m tz^ro*) t ne quadrupole 
components of the "intrinsic" CMB fluctuations. As a 
consequence of statistical isotropy, the a\ m can be as- 
sumed to be equal up to a phase, i.e., a\ a — Qj, 

4i = v^Tse 4 " 1 C/i a 2 
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FIG. 1. The quantity |56o (170) — 5ao(??o)|/«i plotted as a func- 
tion of (1. 



The quadrupole amplitude is denned by Q 



SFEm|C + 4| ■ If ' $ and (^o(%) r <J«o(r?o)) 
are appropriately chosen, then one can obtain a "low" 
value of Q2 — 5.3 x 10~ 6 even if the prediction for a 
spherical Universe, Qi ~ 1.3 x 10~ 5 , would be higher. In 
particular, from the analysis of |24j . we can extract the 
required value of |56o(^o) — Sao(vo)\ — 2 x 10~ 5 . 

It is however important to note that the cancellation 
of the intrinsic quadrupole component by the integrated 
Sachs- Wolfe effect associated to the ellipsoidal Universe 
requires a correlation between the primordial and the late 
components, as discussed already in [18] and quantified 
in [25J . Even in the absence of correlations with the 
intrinsic quadrupole, the effect on the CMB provides the 
strongest constraints on our scenario: those constraints 
can be expressed as \Sb (r] ) — 6ao(rjo)\ < 2 x 10~ 5 . In 
figure [lj we show the log-plot of \5bo(rjo) — 5ao(r]o)\ /k\ 
as a function of the parameter p. This plot shows how, 
for m > O(few) x Hq, the degree of anisotropy grows ex- 
ponentially with m. Note also that, even for m = 0, we 
can have a small anisotropy imprinted on the CMB. As 
we will discuss below, however, in this case the original 
gradient in <fi would not have been negligible during infla- 
tion, and its effect should have been taken into account 
in the calculation of the primordial spectrum of density 
fluctuations. 



V. ON THE INITIAL GRADIENT 

In this section we discuss the magnitude of the initial 
gradient ki(0) = K\ H Mp\. Since for H ^> m the evo- 
lution of the field <f> is frozen by Hubble friction, we have 
cj) ~ Ki(0) z inside our Hubble patch. As a consequence, 
the energy p^ = Tf t scales as K\ (0) 2 / (2 ao(t) 2 ) . This gra- 
dient term redshifts as a(t)~ 2 , i.e., as spatial curvature. 



During inflation, with Hubble parameter Hj, 
a fraction 
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of the background inflaton energy, where N is the num- 
ber of e-foldings from the end of inflation and a cn d is 
the scale factor at the end of inflation. In particu- 
lar, if inflation lasted A; e-foldings and we denote by 
Trh — (3 Afpj i/f) 1 / 4 the reheating temperature (we as- 
sume instantaneous reheating and ignore the effects of 
<7* , the effective number of relativistic degrees of freedom 
in the system), then we can trade «i(0) for (p<f>/ ptot)\i 
computed at the beginning of inflation 




(34) 



where we have used a en d = Tq/Trh with T ~ 2x 10 _4 eV 
being the current value of the CMB temperature. 

As discussed in section IV, data require \Sbo(rj — 
<5oo( ? 7o)| <2x 10 -5 , that implies 



Ptot 



< 10" 



|^o(?/o) - 5ao(vo)\ 



it2 n,f2 
J ± RH 



(35) 

where the quantity \Sbo(rjo) — 5ao(rjo)\ / k\ is plotted in 
figure [T] 

Finally, we use the relation e N ° hs ~ (if/ T )/(H T RH ) 
(iV b s is the number of observable e-foldings of inflation) 
to obtain the condition 



3 x 10" 



\Sbo(f] ) - Sa (r/ )\ 



c 2(W i -7V obs ) 



Ptot 

Therefore inflation can start with a sizable (but still 
smaller than unity) value of (/?^/ptot)|i provided A; is 
sufficiently larger than A b s - Then when the observ- 
able scales left the horizon, A D b s e-foldings before the 
end of inflation, the perturbation in the scalar <j> was 
a smaller fraction (by a factor exp [2 (A Q b s — A^)]) of the 
background energy, small enough to affect neither the dy- 
namics of the inflaton nor that of the perturbations. The 
gradient in <fi stays irrelevant until the Hubble parameter 
is comparable with m, when it will start increasing un- 
der the effect of the tachyonic potential and will give a 
perturbation of the right amplitude today. 

To fix ideas we can set A Q b s = 60, fi = 10 (so 
that \Sb (r)o) - Sao(r)o)\ /«? - 10 5 ), A ; = 71. Then, 
(/V/Ptot)|i — 0.1. At the time observable scales exited 
the horizon, p^j Ptot would have redshifted by a factor 
e 2x (60-71) ^ iq-w an( j ^ wou i c j have been negligible. 
In the absence of the tachyonic mass, the gradient of 
<j> would be still negligible today. However, due to the 
tachyonic enhancement |<56o(??o) — <^o(?7o)|) the effects of 
such a primordial gradient will manifest themselves in 
the lowest CMB multipoles [27] . 

Before concluding this section, let us discuss the pos- 
sible worry that, since <f> is effectively a massless field 
for most of the history of the Universe, its quantum 
fluctuations generated during inflation may be larger 
than the classical gradient that is central to our anal- 
ysis. It is easy to choose the inflationary parameters 
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in such a way that the quantum fluctuations are sub- 
dominant with respect to the original gradient. In- 
deed, the amplitude of the quantum fluctuations in </> 
will be of the order of Hj, the Hubble parameter dur- 
ing inflation. During inflation, the modulation of <f> at 
scales z ~ Hq 1 , which are relevant today, is of the or- 
der of ki(0)2! =s SiMpi ~ v^/^fe^-^Mp!. 
Assuming (p<f> / ptot)\i — 0(1), we see that, as long as 
Hi < Mpje^ *"-^' (e.g., Hi < M P1 in the ex- 

ample ~N\ — 7V bs = 11 considered above), the effects of 
quantum fluctuations are negligible. 



VI. CONCLUSIONS AND DISCUSSION 

If inflation lasted a relatively short time, then some 
primordial gradients could exist as a relic of the chaotic 
pre-inflationary dynamics. A gradient in the quintessence 
field, even if too small to leave any detectable effects 
during the observable epoch of inflation, might be am- 
plified by a tachyonic quintessence to have observable 
magnitude today. We have seen that in the simplest sce- 
nario, the dominant effect of the amplified gradient is on 
the CMB quadrupole, whose observed amplitude sets the 
strongest constraints on the parameters of the problem. 

The fact that CMB fluctuations are affected only 



at the quadrupole level is a consequence of the ap- 
proximation <j) ~ Kx(t) z. Higher powers of z in the 
expansion of (f> would lead to the generation of higher 
multipole contribution. Generically, terms of order z n 
in the expansion of the scalar field <f>{t, z) will generate 
contributions on both Sa(t, z) and Sb(t, z) up to order 
z n_1 , which, in turn, provide terms up to order z n+1 in 
Cisw as for equation (27 1. That is, terms of order z n in 
<p[t,z) will affect Ci's of £ = n + 1. Given the hints of 
multipole alignments up to i ~ 40, it would be interest- 
ing to be able to extend our mechanism in such a way 
that larger powers of z in the expression of 4>{t, z) are 
generated. One possibility is that the self-interactions 
of <f> during the recent cosmological evolution are re- 
sponsible. For instance, if the quintessence field <fi is 
given by a pseudo-Nambu-Goldstone boson, its potential 
is 1/(0) = A [1+ cos (<£//)] /2 (/ 2 = A/(2m 2 )). By 
following our analysis of section |Hj we see that, because 
of the Klein-Gordon equation, the coefficients of lower 
powers in z will act as sources for those of higher 
powers, hence generating higher multipoles. A further 
investigation of this mechanism is subject of future work. 
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